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Abstract. Detailed measurements of the azimuthal anisotropy (2) for identi- 
fied charged particles are reported as a function of transverse momentum (jpt) 
. and centrality for Au+Au collisions at ^/snn =200 GeV. The measurements 

indicate clear evidence for eccentricity and particle flavor scaling over a broad 
. range of centralities and transverse rapidity yx, indicating a hydrodynamical 

origin of the fine structure of azimuthal anisotropy at RHIC. The observed 
scaling supports the picture of a suddenly hadronizing (recombining) fluid of 
fj"^ ' quarks. An apparent breaking of flavor scaling at relatively large values of yj> 

""^ . points to an important change in the mechanism for particle emission. 
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1. Introduction 

X 

During the early stages of an ultra-relativistic heavy-ion collision, an extremely 
high energy density system, possibly consisting of a new phase of nuclear matter, 
is expected to be formed [^ El El EJ El- The dynamical evolution of this matter is 
predicted to reflect its properties El El- Consequently, much effort is currently 
centered on the study of reaction dynamics at the Relativistic Heavy Ion Collider 
(RHIC). Azimuthal correlation measurements constitute an important probe for 
reaction dynamics. They serve as a "barometric sensor" for pressure gradients 
developed in the collision and hence yield insight into crucial issues of thermalization 
and the equation of state (EOS) [El E|- They provide important constraints for 
the density of the medium and the effective energy loss of partons which traverse 
it [ El I10| . They can even provide valuable information on the gluon saturation 
scale in the nucleus [ITT|. 

Azimuthal correlation measurements show significant harmonic strength at mid- 
rapidity with characteristic dependencies on pt and centrality [ 1121 1131 ^] ^1- The 
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anisotropy of this harmonic pattern is typically characterized by the second order 
Fourier coefficient, 



v 2 = (e^* 1 -** J ' ) ), (1) 

where 4>i represents the azimuthal emission angle of an emitted particle and $>rp 
is the azimuth of the reaction plane. The brackets denote statistical averaging over 
particles and events. A large amount of available data for Au+Au collisions at 
,/sWn = 62.4, 130, and 200 GeV [El EI El El E! indicate that the magnitude 
and trends of 2 (for py<2.0 GeV/c) are under-predicted by hadronic cascade models 
supplemented with string dynamics [I18|. but are well reproduced by models which 
incorporate hydrodynamic flow [H Ej. This has been interpreted as evidence for the 
production of a thermalized state of partonic matter 121 El- If this is indeed the 
case, then the fine structure of azimuthal anisotropy (ie. its detailed dependence 
on centrality, transverse momentum, particle type, higher harmonics, etc) should 
reflect the scaling "laws" predicted by ideal hydrodynamics. 

In this work we use detailed differential 2 measurements to test for such scaling 
"laws" and the onset of competing mechanisms. 

2. Hydrodynamic scaling 

An important scaling prediction of hydrodynamic theory is exemplified by the exact 
analytic hydro solutions [El exploited in the Buda-Lund model [[2| . For harmonic 
flow, the model gives: 

V2n = JJ-T> n = M' w = j=r (jr - Tfr) , (2) 
Io{w) 4m t \T y T X J 

where 7o,n are modified Bessel-functions, m t is an average of the rapidity dependent 
transverse mass (at mid-rapidity, Tot = tut), and T x and T y are direction (x and 
y) dependent slope parameters: 

T x = T Q +m t xj n _ , (3) 
1 T + ra t a 

T y = T + m t Y? . (4) 

1 T + m t a 

Here, Xf and Yf gives the transverse expansion rate of the fireball at freeze-out, 
and a = (To — T s )/T s is its transverse temperature inhomogeneity, characterized by 
the temperature at its center To, and at its surface T s . The important prediction of 
Eq. [3 is that the relatively complicated fine structure of azimuthal anisotropy can 
be scaled to a single function. An illustration of this fact can be made for particle 
flavor scaling via substitution of the transverse rapidity [I21j. yx = sinh~ x (j>t j m) 
in Eq. EJto give 

2 m / k 2 T k 3 T 2 \ 

v 2 = kiy T — 1 + + - ...J, (5) 

Jo V Kim k\ m ) 
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where, kl, k2, k3... are largely governed by the expansion rate. Close inspection of 
the leading term in Eq.[S] indicates that V2 for different particle species should scale 
with yfj? — k m • y\m. Here k m is a mass dependent factor which is ~ 1 for relatively 
heavy particles. It is straightforward to show that hydrodynamics also predicts that 
V2 should scale with the spatial eccentricity e — (Y 2 ~ X 2 ) / (Y 2 + X 2 ) , of the overlap 
between the two colliding gold nuclei and the higher harmonic V4 ~ \v\ + k m y\. 

3. Analysis 

The present analysis is based on ~ 22 M minimum bias Au+Au events obtained 
with the PHENIX detector at ^snn — 200 GeV during the second running period 
(2002) at RHIC. Charged tracks were detected in the east and west central arms 
of PHENIX [El EE], each of whicn subtends 90° in azimuth 0, and ±0.35 units of 
pseudo-rapidity r\. Track reconstruction was accomplished at each collision energy 
via pattern recognition using a drift chamber (DC) followed by two layers of multi- 
wire proportional chambers with pad readout (PCI and PC3) located at radii of 
2 m, 2.5 m and 5 m respectively [El EST The collision vertex z along the beam 
direction was constrained to be within \z\ < 30 cm. A confirmation hit within a 2a 
matching window was required in PC3 and the electromagnetic calorimeter (EMC 
PbSc) or the time-of-flight detector (TOF), to eliminate most albedo, conversions, 
and decays. Particle momenta were measured with a resolution of Sp/p — 0.7% © 
1.0% p (GeV/c). 

Event centralities were obtained via a series of cuts in the space of BBC versus 
ZDC analog response; they reflect percentile cuts on the total interaction cross 
section [ESj- Estimates for the number of participant nucleons N part , were also 
made for each of these cuts following the Glauber-based model detailed in Ref. [ 

E3|. 

In this analysis the combination of the TOF detector and six sectors of the 
(EMC PbSc) was used to identify charged particles. Particle time-of-flight was 
measured using the TOF ( or EMC ) and the collision time defined by beam counters 
(BBC). A timing resolution of ~ 120 ps and ~ 370-400 ps was obtained for the TOF 
and the EMC (PbSc) respectively. This allowed meson ( 7r ± , K ± ) and baryon ( p , 
p ) separation up to a pr — 3.5 GeV/c with the TOF detector and up to 2.5 GeV/c 
with the (EMC PbSc). 

The differential «2(pt, centrality) measurements were obtained via the reaction 
plane technique which correlates the azimuthal angles of charged tracks detected 
in the central arms with the azimuth of an estimated event plane $2, determined 
via hits in the North and South BBC's located at | 77 |~ 3 - 3.9 [E33 E3- Due the 
large rapidity gap (~ 3 — 3.9 units) between the particles used for reaction plane 
determination and the mid-rapidity particles correlated with this plane, one expects 
that the latter correlations are less influenced by non-flow contributions especially 
for pt < 3.0 GeV/c. The estimated resolution of the combined reaction plane 
from both BBC's [ 1151 EH has an average of 0.33 over centrality with a maximum 
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of about 0.42 for Au+Au collisions at ^/snn~ 200 GcV. Thus, the estimated 
correction factor, which is the inverse of the resolution for the combined reaction 
plane, ranges from 2.4 to 5.0. 

3.1. V2 results 

Figures ^ and [2] summarize the centrality and pt dependence of «2 for charged 
mesons( tt^, K ± ) and baryons ( p + p ) detected in the EMC+TOF respectively. 
They give an excellent overview of the the evolution of i>2 as centrality and px 
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Fig. 1. V2 vs pt for charged mesons ( n^, K ) for different bins in reaction 
centrality (top to bottom) 0-5%, 5-10%, ... , 35-40%, 40-50%, 50-60%. The error 
bars shown indicate statistical errors only. Systematic errors are estimated to be 
less than 10%. 



are varied. The results shown for protons and anti-protons give an especially good 
view of the evolution away from the well know quadratic dependence of V2(pr) 
(which is also observed in very central collisions for these data) as the collisions 
become more peripheral. Such a dependence could result from changes in the freeze- 
out temperatures and/or the radial flow velocity as the collisions become more 
peripheral. 
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Fig. 2. V2 vs pt for charged baryons ( p + p ) for different bins in reaction 
centrality (top to bottom) 0-5%, 5-10%, ... , 35-40%, 40-50%, 50-60%. The error 
bars shown indicate statistical errors only. Systematic errors are estimated to be 
less than 10%. 



3.2. Eccentricity scaling 

As indicated earlier, hydrodynamics [ \2'6\ predicts eccentricity scaling of the az- 
imuthal anisotropy. To test for this scaling, one can divide the vi values obtained 
at a given centrality by the eccentricity e obtained from a Glauber-based calcula- 
tion obtained for the same centrality Alternatively, one can simply divide the 
differential anisotropy V2(N par t,pr) obtained at a given centrality by the pr inte- 
grated value for the same centrality selection V2{N part ). The underlying idea here 
is that the integral flow is monotonic and linearly related to the eccentricity over 
a broad range of centralities [EH- Another advantage of this approach is that the 
ratio V2(N part , pr) /v2{N part ) gives a scale invariant variable which automatically 
reduces the systematic errors associated with an eccentricity evaluation, and the 
reaction plane resolution. 

Fig. |21 shows the V2 of charged mesons obtained for different centralities scaled 
by the value of the integral flow obtained for each of these centralities. The figure 
shows essentially perfect scaling for mesons as would be expected for a process driven 
largely by the eccentricity of the overlap region of the two colliding nuclei. The 
results of a similar scaling test for baryons are shown in Fig.^J Despite a continuous 
evolution in the shape of V2 vs pr with centrality, they indicate a relatively good 
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Fig. 3. V2 vs pr for charged mesons ( for different bins in reaction centrality ) 
scaled by the value of integral flow for each centrality. . 



scaling for the two centrality ranges: 0-20% ( left panel ) and 20-50% ( right panel ). 



3.3. Flavor scaling 

Following Eq.0 the scaling variable = k m ■ y^m, was used to investigate flavor 
scaling. The results of such a test is summarized in Fig. [5] The left panel of the figure 
shows a comparison of V2{pt) for protons, kaons and pions measured with the TOF 
detector for the centrality selection 5-30%. The characteristic flavor dependence of 
vi is cleanly evidenced by these data, ie. mass ordering at low momentum and a 
reversal of the magnitudes of the v% values for baryons and mesons for pr ~ 1 .8-2.0 
GeV/c. The right panel of Fig. [3] and Fig. show that very good scaling of v-i is 
achieved with yt^ — k m ■ y^m in accordance with the predictions of hydrodynamics. 
Fig. gives an illustration of y£ s scaling for combined results which include the 
neutral kaons and lambda hyperons measured by the STAR collaboration. The 
latter were obtained in Au+Au collisions at ^snn =200 GeV for the same centrality 
selection [El- Although the data shown in Fig.HJindicate rather good scaling for all 
measured particles over a broad range of y£ s , one can see clear evidence for a break 
in this scaling for y^f > 2. Such a break could be signaling a change in mechanism 
for high pt particles or a break down of ideal hydrodynamic flow. 
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Fig. 4. V2 vs pt for charged baryons ( for different bins in reaction centrality ) 
scaled by the value of integral flow for each centrality. Left panel corresponds to 
centrality range 0-20% and right panel to centrality range 20-50% 
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Fig. 5. V2 vs pt for charged pions, kaons and protons (left panel). The right panel 

shows the result of the y T s scaling predicted by hydrodynamics. Results are shown 
for the centrality selection 5-30%. 



3.4- Higher harmonic scaling 

As was pointed out earlier, ideal hydrodynamics [ [20! gives predicts a straightfor- 
ward scaling relationship between the second harmonic i>2 and the higher harmonics 
of azimuthal anisotropy ( uj, vq, ...). It is easy to show that the leading term of the 
relationship between vi and can can be expressed as 

V4 = ^ • V 2 + k 5VT> ( 6 ) 

The measurement of the higher harmonics of azimuthal anisotropy are not available 
in PHENIX ( ongoing ). However, we can test the predicted scaling relationship 
between v-i and via data published by the STAR collaboration [|2H1 129j . Figure 
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Fig. 6. y£ s scaling of charged pions , kaons and protons (PHENIX) and neutral 
kaons and lambdas (STAR) [Ej- Results are shown for the centrality selection 
5-30%. 



[7| shows the (y^) 2 dependence of 1)4 for charged particles from Au+Au collisions 
at y/SNN =200 GeV and V2 values scaled according to Eq. The figure indicates 
rather good agreement between the measured V4 values and scaled V2 values over a 
broad range in ■ 

4. Summary and Conclusions 

In summary detailed measurements of the fine structure of elliptic flow in Au+Au 
collisions at y/s^N =200 GeV measured by PHENIX collaboration at RHIC are 
presented. They show eccentricity scaling and flavor scaling over a broad range 
of centralities and particle flavors. This observed scaling gives strong support for 
essentially ideal hydrodynamical flow at RHIC. The observed scaling also supports 
the picture of a suddenly hadronizing (recombining) fluid of quarks. 
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